Optimal cooling of a driven artificial atom in dissipative environment 
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We study microwave-driven cooling in a superconducting flux qubit subjected to environment noises. For 
the weak decoherence, our analytical results agree well with the experimental observations near the degeneracy 
point and show that the microwave amplitude for optimal cooling should depend linearly on the dc flux detuning. 
With the decoherence increasing, more vibrational degrees of freedom couple in, making the ordinary cooling 
method less effective or even fail. We propose an improved cooling method, which can eliminate the perturba- 
tion of additional vibrational degrees of freedom hence keep high efficiency even under the strong decoherence. 
Furthermore, we point out that the decoherence will modulate the frequency where microwave-driven Landau- 
Zener transition reaches maximum in both methods, displaying the feature of incoherent dynamics which is 
important for the optimal cooling of qubits and other quantum systems. 

PACS numbers: 37.10.De,03.65.Yz,85.25.-j,03.67.Lx 



I. INTRODUCTION 

Superconducting devises based on Josephson tunnel junc- 
tions, shown to act as artificial atoms, can be used to demon- 
strate quantum phenomena at macroscopic scales and hold 
promise for applications in quantum computation as qubits 
Qubits are very delicate and their manipulation with elec- 
tronics needs to be isolated from the disturbance of the en- 
vironmental thermal noisei^"— Therefore, cooling the qubits, 
considered as a straightforward means to preserve the coher- 
ence of qubits, paves the way for remarkable achievements in 
demonstrating various quantum coherent phenomena. In gen- 
eral, the dilution refrigeration is used for cooling with some 
drawbacks such as limited cooling efficiency and poor heat 
conduction. Therefore, novel cooling techniques are required 
to improve the performance of qubits. 

Recently an interesting technique, microwave-driven cool- 
ing, was developed from the extensive investigations on cou- 
pled systems of the qubits and quantum resonators. It was 
found that the qubit can cool down the resonator—"—. On the 
other hand, in a recent experiment, the inverse process, us- 
ing a "resonator" to cool a qubit, was also proposed to lower 
the temperature of a superconducting flux qubit^^ by Valen- 
zuela et al}^ . At optimal cooling, they achieved the effec- 
tive temperatures 3 mK, which was two orders of magnitude 
lower than the surrounding environment temperature, with the 
condition of the amplitude where the particular flux ampli- 
tude made only one of two side crossovers reached and the 
frequency of 5 MHz. This process is analogous to the op- 
tical sideband cooling,^^"^^ where a two-level system (TLS) 
in an ion (an "internal" qubit) is used to cool ion trap po- 
tential (an "external" oscillator), subsequently providing di- 
rect cooling to the degrees of freedom of interest. But dif- 
ferent from the optical method, here an ancillary "internal" 
oscillator-like state is coupled to an "internal" qubit by a side 
crossover, cooling the "internal" qubit through Microwave- 
driven Landau-Zener (MDLZ) transitions (or Landau-Zener- 
Stiickelberg interferometry)J^i^2i^ As the qubit is driven 
through a degeneracy point (or crossover), Landau-Zener (LZ) 
transition happens which usually is used to enhance the 
quantum tunneling rater^ to prepare the quantum statesr^ 



and to control the gate operations. ~ When a periodic strong 
driving is applied to the qubit, the qubit undergoes repeated 
LZ transitions at the degeneracy point. If the driving fre- 
quency is larger than the decoherence rate, the repetition of 
LZ transitions can induce quantum-mechanical interference, 
which leads to an oscillatory dependence of qubit popula- 
tion in the final state on the detuning from the degeneracy 
point and the microwave amplitude, known as Landau-Zener- 
Stiickelberg (LSZ) interferometry^-i2i21^-i^ When the deco- 
herence rate is larger than the driving frequency, the interfer- 
ence patterns are nearly washed out. Due to the loss of coher- 
ence, this incoherent region draws little attention in quantum 
information until the microwave-driven cooling emerges^ 

In this paper, firstly we obtain the optimal cooling con- 
ditions in Ref. [27] based on a microscopic model of LZS 
interferometry in the multi-level system.^^ Besides only one 
crossover reaches, the microwave amplitude for optimal 
cooling should be linearly dependent on the dc flux detuning. 
Moreover, the frequency condition also varies slightly with 
different detunings. Then we investigate how the decoherence 
modifies the condition for optimal cooling. The strong deco- 
herence might couple the "internal" qubit to another "inter- 
nal" oscillator-like state which would interact with the origi- 
nal oscillator-like state, destroying the cooling effect although 
controlled amplitude makes only one of two side crossovers 
reached. To realize the cooling under the strong decoherence, 
we propose an exquisite microwave manipulation to exclude 
the excessive oscillator-like state which demolishes the cool- 
ing in MDLZ transitions. Active cooling of qubits provides 
an effective means for qubit state preparation even under the 
strong decoherence and for suppressing the decoherence in 
multi-qubit systems. Especially, it can be applied to refrig- 
erate environmental degrees of freedom"^^ and cool neighbor- 
ing quantum systems. Additionally, for systems consist of 
the interaction between the qubit and the real resonator, such 
as micromechanical beam, this method can also be used to 
cool down the qubit, or even resonators. We also find that 
the frequency for optimal cooling relies on the decoherence, 
not always achieved at the extremely low frequency. The de- 
coherence will modulate the frequency corresponding to the 
maximum rate of MDLZ transition, leading the change of the 
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FIG. 1: (Color online). Schematic energy diagram of a multi-level 
superconducting flux qubit. (a). Red solid curve represents the 
microwave $ac in the ordinary cooling method in Sec. III. Black 
dashed curve represents the microwave in the new cooling method 
in Sec. IV. The dashed black line expresses a particular static flux 
detuning 6^dc- The symmetry makes |$i2| = l^osl- States |1) 
and 1 3) are in the left well; |0) and |2) are in the right well. The 
red solid path describes the cooling with the population in state |1) 
transferred to |0). The blue dot path represents the population trans- 
ferred from state |0) to (b). The shaded region describes the 
"internal" qubit. Conventionally, we employ a microwave to couple 
a oscillator- like state |2) to the "internal" qubit. However, the strong 
decoherence would make another oscillator- like state |3) destroys the 
cooling completely, although state |3) has little effect under the weak 
decoherence. 



frequency condition in both methods. At the strong deco- 
herence environment, LZS interference displays the station- 
ary population of incoherent evolution. While the trend of 
the quantum information is coherent manipulation, our analy- 
sis shows the importance of the incoherent LZS dynamics in 
quantum cooling. 

The organization of this paper is as follows. In Sec. II, 
we give the basic model of LZS interferometry, including LZ 
transition and the rate equation. Then we divide the cooling 
region into two regions and the incoherent transition region is 
introduced. In Sec. Ill, we present a detailed analysis of the 
optimal cooling in the experiment and find that there are dif- 
ferent amplitude conditions for optimal cooling at different 
detunings, besides the change of the frequency condition. In 
Sec. IV, we study the cooling under the strong decoherence, 
and find that the strong decoherence makes the cooling inef- 
fective in the ordinary method. Then we propose an improved 
cooling method to solve the problem and discuss the optimal 
cooling in the new method. In Sev. V, we study the relation 
between the efficiency of MDLZ transition and the decoher- 
ence, explaining abnormal phenomena. 



II. BASIC MODEL 

We focus on a multi-level superconducting flux qubit, 
which consists of a superconducting loop interrupted by three 
Josephson junctionsi2i22ill (Fig. 1). If the external flux bias 
$ ^ 0.5<l>o, where = Z?^/ 2 e is the flux quantum, a double 
well landscape of the potential energy parameterized by the 
dc flux detuning S^dc = ^ — O.b^o exists in the flux qubit. 
In the millikelvin temperature ~ 10 mK, a series of diabatic 
states, denoted as \i) and | j), respectively are localized in dif- 
ferent wells with the crossovers Aij between states \i) and \j) 
(i = 0, 2, right well with negative slope; j = 1, 3, left well with 
positive slope). We assume the barrier is relatively high thus 
the tunneling strength (A^j) is small. The two wells corre- 
spond to different flowing directions of the persistent currents, 
namely clockwise and anticlockwise. S^ij is the flux detun- 
ing, at which the crossover Aij is reached. As the qubit is 
driven with a microwave ^ac = ^rf sincjt, where ^rf is the 
microwave amplitude, the time dependent flux detuning is 

5^{t) = S^dc + ^ac = ^^dc + ^rf smut. 

Then the time dependent energy of eigenstate \i) can be de- 
scribed as 

e^{t)=e;±m^S^t) 
= el ± niiib^dc + ^r/ sincjt), 

where e\ is the energy spacing between state \i) and the lowest 
state in the same weU at b^(t) = m^o, £i = ± mi6^dc 
is the dc energy of eigenstate |^), (+ corresponds the states in 
the left well, — corresponds the states in the right well), and 
nii = dEi{^)/d^ is the energy level slope of 

When the state |z) is driven through the crossover A^j at the 
time t = to, the probability of LZ transition from state \i) 
\j) is given by (we set h = ks = I) 

A2. 

p = l-exp(-^), 
where the relative-energy sweep rate ( = ^[\ei{t) 



to 



(1) 



t = to- 



For the case of the periodic sinusoidal driv- 



ing, we have \€i{t) — ej{t)\ = {\nii\ + |mj |)<l>^/ sincjt and 
C = (l^il + \mj\)uj^rf\cosujto\. Therefore C oc u, which 
means larger frequency leads to less probability of LZ transi- 
tion. 

To describe a driven multi-level qubit system interacting 
with a dissipative environment, we consider the Hamiltonian 



H — Hq + 



int 



(2) 



where Hq = %^K)0'l ^T.^k{t)\k){k\ is the Hamilto- 

i^j k 

nian of the multi-level qubit. Hint is the system-environment 
interaction Hamiltonian, and Hb is the environment Hamilto- 
nian. Considering the noise from the environment, we assume 
the linear coupling between the qubit operators and environ- 
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ment operatori^iiS 

Hint = -f E (10(^1 -U)(J|) + 
|i)Gright well 
|j)Gleftwell 

Y,^kk'(\k){k'\^W){k\)Q, (3) 

where Q is the operator of environment noise, states \k) and 
\k') belong to the same well and (j)kk' is the coupling matrix 
element. The first term describes the dephasing from low fre- 
quency noise. The second term Hhf = J] (j)kk' {\k){k'\ + 

k^k' 

\k'){k\)Q represents the effect of high frequency noise to pro- 
duce the intrawell relaxation. 

To describe the long-time system evolution, we can employ 
the rate equation approach. In this case, MDLZ transition 
between two states in opposite wells needs to be considered. 
Starting from the Hamiltonian (2), we use the Gaussian white 
noise model and perturbation theory^^i^i^ to derive MDLZ 
transition rate from state |z) to | j) as 

ly. . - ^ V (r2+7)j;(A,,/u;) ... 
n 

where Sij = Si — Sj is the dc energy detuning of states \j) 
and Aij = \mi\^rf + l^j l^r/ is the energy amplitude of 
the microwave, r2 is the dephasing rate, 7 = {Tji + r^o)/2, 
Tji and r^o are the intrawell relaxation rate from states \j) 
and \i) in the left and right well respectively. Of course, when 
i=0 or j=l, the corresponding F^o or F^i is zero. The term 
r2 +7 describes the decoherence in this sub-system composed 
of states \j) and In this article, when considering differ- 
ent decoherence, we would hold the intrawell relaxation rate 
Tai = = 27r X 0.1 GHzf^ only changing and mentioning 
the dephasing rate. 

In order to analyze the optimal cooling, we consider four 
lowest states in which the qubit state occupations pi{i = 
0,1, 2, 3) follow 

Poo = -Pooi^oi + 1^03 + Wol) + (r^o + Wio)pii^ 

^30P33 + r20P22 

Pii = -Pii(r'io + W12 + Wio) + (Foi + 1^01)^00 

+r3iP33 + W21P22 
P22 = -P22(r20 + W21 + W23) + "^^32^33 + 1^12^11 

Poo + Pll + P22 + P33 = 1, (5) 

where F^q, Fq^ = F^q exp(— £oi/^) is the down and up in- 
terwell relaxation rate between states |1) and |0), and T is the 
temperature of the environment. 

In the microwave-driven sideband cooling, states |0) and 
|1) are the two lowest states, which are far below the other 
states and form an "internal" qubit. The state |2) is treated as 
an ancillary "internal" oscillator-like state, which is localized 
in the same well with state |0), and A 12 is considered as the 
side crossover. If the qubit is in equilibrium with the envi- 
ronment, the population in state |1) is excited from state |0) 
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FIG. 2: (Color online). Schematic energy diagram illustrating the 
cooling regions. In (a), the solid line expresses the transition from 
state 1 1) to 1 0) at the coherent region. The arrows mark well-resolved 
resonances of n-photon transition observed under the high frequency 
driving. The shaded region expresses the interference phase. In (b), 
the coherence interference vanishes and the transition from state |1) 
to |0) is realized at the incoherent transition region. 



thermally following the Boltzmann relation 

Pii/poo = ex.p{-eio/T). (6) 

This expression shows the influence of temperature is large 
at small detunings near m<l>o, which thereby becomes main 
cooling interests. As shown in Fig. 1(a), by repeated LZ tran- 
sitions, the microwave drives the unwanted thermal popula- 
tion in state 1 1) to |2), whose population relaxes fast into state 
|0). Then the population transferred to state |0) faster than the 
thermal repopulation of state |1), cools the "internal" qubit. 
As indicated from Eq. (6), less population in state |1) corre- 
sponds to lower effective temperature. Hence we can use the 
population in state |1) to determine the effective temperature 
realized by microwave cooling. 

Microwave-driven sideband cooling depends on MDLZ 
transition from state |1) to |2) which exhibits a rich structure 
as a function of the microwave frequency. Therefore we divide 
the cooling into two regions by the frequency. As shown in 
Fig. 2 (a), at the coherent region where the frequency is larger 
than the decoherence rate F2 +7, the qubit state is split into the 
coherent superposition of states |1) and |2) by the crossover 
A 12. Then a phase difference is constructed to form the quan- 
tum interference.^^ Especially at high frequencies near several 
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FIG. 3: (Color online). The population in state |1) versus flux de- 
tuning and microwave frequency with the flux amplitude $r/ = 
8.4 m$o- The locations of the crossovers Aoi,and A12 are 
m$o, and 8.4 m$o, respectively. The diabatic energy-level slope 
|mo|(|mi|) =27r x 1.44 GHz/m$o, and |m2|(|m3|) =27r x 1.09 
GHz/m$o. Aoi/27r = 0.013 GHz, Ao3/27r = Ai2/27r = 0.09 
GHz. The parameters above are from the experiment.— Further- 
more we estimate A23/27r — 0.5 GHz. The dephasing rate we 
used in calculation r2/27r = 0.06 GHz, the intrawell relaxation 
rate r2o/27r = r3i/27r = 0.1 GHz,— and the interwell relaxation 
rio/27r = 0.00005 GHz. The environment temperature T is 50 mK. 
(b). The population in state |1) versus the microwave frequency at 
the detuning 0.05 m$o- (c) The population in state |1) versus flux 
detuning at the frequency 1460 MHz. 



subgigahertzes, the well resolved n-photon resonances are ob- 
served. At the incoherent region with the frequency is lower 
than the decoherence rate, as state |1) is driven through the 
crossover A 12, the state is split into the coherent superposi- 
tion of states |1) and |2) by LZ transition. However the deco- 
herence eliminates its coherence quickly such that the super- 
position becomes the incoherent mixture, as illustrated in Fig. 
2 (b). When the microwave drives the state back to A 12, the 
second LZ transition occurs, but state 1 1) split from state 1 1) 
would not interfere with state 1 1) split from state |2) due to the 
lost coherence. In next sections we will show this region has 
a important role in the optimal microwave cooling. 



III. OPTIMAL COOLING UNDER THE WEAK 
DECOHERENCE 

In the experiment under the weak decoherence,^^ e. g. 
r2/27r = 0.06 GHz and 7/27r = 0.05 GHz, Valenzuela et al 
proposed that the optimal cooling near the degeneracy point 
was realized at the flux amplitude ^rf = ^12 = 8-4 m<l>o and 
the frequency ~ 5 MHz. 
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FIG. 4: (Color online), (a) The population in state |1) versus mi- 
crowave frequency and flux amplitude with the detuning 6^dc — 
0.05 m$o- The other parameters used are the same with those of 
Fig. 3. (b) The lowest population extracted from (a) at each fre- 
quency versus microwave frequency, (c) The amplitude correspond- 
ing to the lowest population extracted from (a) at each frequency 
versus microwave frequency. 



In the stationary case, P22 = Pii = P33 = Poo = 0. with 
the amplitude ^rf =8.4 m^o we use the rate equation Eq. 
(5) to plot the population in state |1) at different detunings 
and frequencies as shown in Fig. 3 (a). At the coherent region 
the effective cooling can only be achieved at some detunings, 
[Fig. 3 (b)] such that we do not consider the coherent region in 
microwave cooling. At the incoherent transition region, lower 
frequency leads to less population in state 1 1) at all detunings, 
which means more efficient cooling. These results agree with 
the experimental observations.^^ Then we extract the data in 
the detuning d^dc = 0.05 m<l>o [the black dashed line in Fig. 
3(a)] for detailed analysis in Fig. 3 (c). When the frequency is 
less than 100 MHz, the behavior of population dependent on 
the frequency follows the above discussion. However, as the 
frequency reaches 10 MHz, the optimal cooling at this detun- 
ing is realized (point II). Then at point I where the frequency 
is 5 MHz, the population is larger than that at point II. This 
result suggests that the condition for optimal cooling obtained 
by the experiment is not sufficiently optimized. 

To explain the behavior exhibited in Fig. 3, we make a com- 
prehensive analysis by varying microwave amplitudes and fre- 
quencies at the detuning 5^dc = 0.05 m^o, and show the 
population in state |1) in Fig. 4 (a). In detail, we extract the 
lowest population and the corresponding amplitudes at differ- 
ent frequencies, as shown in Fig. 4 (b) and Fig. 4 (c) respec- 
tively. Seen from Fig. 4 (b), the population in state |1) dis- 
plays a monotonic decrease with the frequency and the lowest 
population at the detuning 5^dc = 0.05 m<l>o is realized at 
point III with 5 MHz and the corresponding amplitude 8.35 
m<l>o, instead of I, or II [Fig. 4 (c)]. As a result, at the de- 
tuning S^dc = O.Ob m<l>o the frequency condition for optimal 
cooling is consistent with that of the experiment, while the 
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FIG. 5: (Color online), (a) The population in state |1) versus flux 
detuning and microwave flux amplitude with the frequency 5 MHz. 
The other parameters are identical with those of Fig. 3. (b) The 
amplitude corresponding to the lowest population extracted from (a) 
at each detuning versus the detuning. 



amplitude condition not. 

Then at 5 MHz, we plot the population versus flux detuning 
and amplitude, shown in Fig. 5 (a). At each detuning, we 
extract the amplitude corresponding to the lowest population, 
as shown in Fig. 5 (b). The effective cooling is realized with 
the amplitude ^rf ^ ^12 - ^^dc, where ^rf + ^^dc ^ 8-4 
m<l>o at the detuning near S^dc = m<l>o and ^rf + ^^dc ^ 
8.45 m^o at other detunings. (quantitive result of cooling can 
be seen in Fig. 8) Therefore, at the detuning approaching to 
m^o, the amplitude condition can be changed to ^rf = *^12, 
just as those demonstrated in Ref. [27]. 

Further, we extend the above discussion to all detunings. 
Firstly, following the similar way to obtain Fig. 4 (b), we 
change the amplitude to obtain the lowest population at each 
detuning and frequency, shown in Fig. 6 (a). At each detun- 
ing, the frequency 5 MHz always produces an efficient cool- 
ing. If we extract the data at the detuning S^^c = 0.5 m<l>o, a 
strange phenomenon that the population at 10 MHz is slightly 
less than that at 5 MHz can be observed [Fig. 6 (b)]. It should 
be mentioned that this case is different from the result in Fig. 
3(c) which occurs at the amplitude 8.4 m^o- The explanation 
would be given in Sec. V. Then we take account of the lowest 
population at other detunings. At these detuning, our calcu- 
lation shows that the frequencies corresponding to the lowest 
population are always near 5 MHz and the difference between 
the population of 5 MHz and the lowest population is small. 
Therefore we can still accept the previous frequency condition 
for optimal cooling in practice. Moreover, the amplitude cor- 
responding to the lowest population at each detuning is shown 
in Fig. 6 (c) and accords with that in Fig. 5 (b). 

In the stationary case, considering W03 <C Fai, r2o, 
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FIG. 6: (Color online), (a) The lowest population in state |1) among 
different flux amplitudes versus the frequency and detuning. The 
other parameters used are the same with those of Fig. 3. (b). The 
population in (a) versus the frequency at the detuning 0.5 m$o- (c). 
The amplitude corresponding to the lowest population extracted from 
(a) at each detuning versus the detuning. 



the population in state 1 1) would be obtained as 



Pii 



- 01 



Wo3 



10 



01 



Wi2 + Wo3 



(7) 



When the dc flux detuning is large, I^i2, Tq^ ^ Wqs such 
that Wo3 can be neglected in Eq. (7) and the population in 
state 1 1) would depend on Firstly, we discuss the ampli- 
tude condition for optimal cooling. MDLZ transition rate W12 
can be treated as the population transferred from state |1) to 
1 2) in one unit time. Therefore W12 depends on the population 
transferred by one time of LZ transition in the crossover A12, 
and the frequency which determines how many times of LZ 
transitions occur in one unit time. It is assumed that in the dc 
flux detuning far away from the crossover A 12 there is nearly 
no tunneling transition between states |1) and |2), which is 
given by transition rate^ 



^12 



(r2 + 7) 



2 s?2 + (r2 + 7)^ 



(8) 



where 7 = r2o/2. Eq. (8) indicates the tunneling can be 
considered in a transition region |£:i2| ^ r2 + 7, where F 
~ A^2/[2(r2 + 7)]. Therefore longer time staying in the 
transition region in a period would contribute to larger popu- 
lation transferred from state |1) to |2) and better cooling effi- 
ciency. The long time staying in the transition region means a 
small velocity. At the same time, the velocity of a sinusoidal 
movement is the smallest when it closes to its full amplitude. 
Therefore, the amplitude which just reaches the transition re- 
gion would bring the largest transition rate from state |1) to 
1 2). When the decoherence is weak, the transition region is 
always near the crossover. As a result, we have the ampli- 
tude relation ^rf ^ *^12 — ^^dc as obtained above. This 
relation also ensures that only the side crossover A 12 reached. 
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FIG. 7: (Color online), (a) The population in state |1) versus mi- 
crowave frequency and flux amplitude with the detuning 6^dc — 
0.05 m$o and the dephasing rate T2/2ti — 1 GHz. The other pa- 
rameters used are the same with those of Fig. 3. (b) The lowest 
population extracted from (a) at each frequency versus microwave 
frequency, (c) The amplitude corresponding to the lowest population 
extracted from (a) at each frequency versus microwave frequency. 



On the other hand, if we increase the frequency, smaller pe- 
riod indicates more times of LZ transitions in one unit time 
so as to cause a larger MDLZ transition rate. However, larger 
frequency results into larger velocity and shorter time in the 
transition region. Therefore it may also lead to smaller pop- 
ulation transferred by LZ transition, shown in Eqs. (1) and 
(8). Consequently, the final result is complex and it is hard to 
determine the frequency condition for optimal cooling which 
would be discussed in Sec. V. 



IV. NEW COOLING METHOD UNDER THE STRONG 
DECOHERENCE 

Although the precise microwave manipulation can avoid 
reaching the another side crossover A03 thus eliminating the 
effect of state |3) in the cooling cycle, the strong decoherence 
may still couple states |0) and |1) to state IS)"^^"^^. State |3) 
would act as another "internal" oscillator-like state and de- 
struct the cooling through MDLZ transition in A03, shown in 
Fig. 1 (b). While the microwave drives the state near the side 
crossover A03, LZ transition would transfer the population 
from state |0) to |3), which then relaxes into state In this 
path, state |3) yields an inverse effect to destroy the cooling. 
The steady population in state |1) would be described by Eq. 
(7). 

Under the strong decoherence, e. g. r2/27r = 1 GHz and 
7 = 0.05 GHz, Fig. 7 (a) is the contour plot of the population 
in state |1) at the detuning 5^dc = 0.05 m<l>o with differ- 
ent amplitudes and frequencies. For each frequency, we ex- 
tract the lowest population and the corresponding amplitude, 
shown in Fig. 7 (b) and (c) respectively. Although even un- 
der the strong decoherence the cooling at this detuning still 
realizes the max efficiency at 5 MHz, the decoherence makes 
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FIG. 8: (Color online) (a) The lowest population in state |1) among 
different flux amplitudes versus the frequency and detuning under the 
dephasing rate r2/27r = 1 GHz. The other parameters used are the 
same with those of Fig. 3. (b). The comparison between the lowest 
population in state |1) at 5 MHz and each detuning under different 
dephasing rates (blue dashed line corresponds to V2/2ti = 0.06 GHz 
while red dotted line corresponds to r2/27r = 1 GHz). The solid line 
represents the population in the equilibrium with the temperature 50 
mK. 



the population near 0.5 and the cooling very inefficient. At 
each detuning and frequency, we change the amplitude to ob- 
tain the lowest population in state |1) as shown in Fig. 8 (a). 
Then we extract the lowest population at 5 MHz from Fig. 8 
(a) and plot them as the red dotted line in Fig. 8 (b). The 
population of optimal cooling under the strong decoherence 
is comparable with the population in the equilibrium, or even 
worse. The blue dashed line represents the lowest population 
at 5 MHz under the weak decoherence with r2/27r = 0.06 
GHz. Therefore, the comparison in Fig. 8 (b) shows that an- 
other oscillator-like state |3) counteracts the effect of original 
oscillator-like state |2). 

At the same time. Fig. 7 (c) shows that the optimal cool- 
ing at the detuning 0.05 m<l>o is realized with the amplitude 
8.1 m<l>o and frequency 5 MHz. In this case, the state does 
not reach the side crossover A 12 and the result violates the 
amplitude condition we obtained in Sec. III. The strong de- 
coherence would result in a large transition region near A 12. 
Subsequently in order to realize better cooling larger ampli- 
tude is needed. However larger amplitude would also drive 
the state closer to A03 . Moreover the strong decoherence en- 
larges the transition region near A03 such that the population 
transferred from state |0) to |3) is not able to be neglected, 
which then relaxes into state |1) and causes increasing popu- 
lation in state Therefore, if the amplitude is decreased to 
minimize the efficiency of state |3), although the microwave 
does not completely cover the transition region near A 12, the 
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FIG. 9: (Color online). The lowest population in state |1) among 
different amplitudes versus dephasing rate and microwave frequency 
at the detuning 6^dc =0.05 m$o- The other parameters used here 
are the same with those of Fig. 3. 



final cooling efficiency may be better. By the way, under the 
weak decoherence, the amplitude condition at the detuning 
near the degeneracy point, e. g. S^dc = 0.05 m<l>o, is differ- 
ent from that at detuning far away from the degeneracy point, 
e. g. ^rf + ^^dc ^ 8-45 m$o in Fig. 5 (b). That is be- 
cause at the detuning near the degeneracy point the condition 
^rf + ^^dc ^ 8-45 m<l>o would also make A03 involved 
hence destructing the optimal cooling. Therefore, the final 
amplitude condition is revised to ^rf + ^^dc ^8.4 m<l>o to 
minimize the effect of Wqs . 

At the detuning S^dc = 0.05 m^o, we change amplitudes 
to extract the lowest population at different frequencies and 
dephasing rates, shown in Fig. 9. Independent of the deco- 
herence, smaller frequency results into better cooling. More- 
over, the dephasing rate larger than 27r x 0. 1 GHz would easily 
destroy the cooling. It should be mentioned that in previous 
studies^i*^ the large amplitude microwave drives the system 
through the crossover A03, coupling the two lowest states to 
state 1 3). However here instead of the large amplitude, the 
strong decoherence couples state |3) to the two lowest states. 
As shown in Figs. (7) and (8), in the case of the strong de- 
coherence, decreasing the amplitude can not keep away from 
the influence of state |3). Therefore, new method is needed 
for the efficient cooling. 

In the ordinary cooling cycle, at the beginning the sinu- 
soidal microwave drives the state in the direction to the side 
crossover A 12 and after going back to the initial detuning, the 
state is also driven in the direction to another side crossover 
A03. Just the latter process brings the coupling to state |3). 
Therefore if we cancel the latter process, the cooling would 
come back. That means when the state goes back to the ini- 
tial position, the microwave drives it still in the direction to 
A 12, not A03. This seems a minor change but brings great 
efficiency to the cooling. It is worth to point out that the tech- 
nique is readily available"^^, making our method practical. We 
also can use periodic triangular waveform to drive the state. 
But to calculate conveniently, we still use the harmonic wave- 
form shown in Fig. 1 (a). Then in Eq. (4), MDLZ transition 
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FIG. 10: (Color online), (a) The population in state |1) versus flux 
amplitude and detuning with the frequency 5MHz and the dephasing 
rate r2/27r = 1 GHz in the original cooling method. The other pa- 
rameters are identical with those of Fig. 3. The white dashed line in- 
dicates the optimal cooling region in the original method. The green 
solid line expresses the optimal cooling region in the new method, 
(b). The population in state |1) versus the flux amplitude and de- 
tuning in the same condition with (a) under the new method, (c) 
and (d). The lowest population and the corresponding amplitude ex- 
tracted from (b) at each detuning versus the flux detuning. 



rate can be changed to 

A 2 

W -W - \^ (r2+7)^;(A.,/u;) .gx 

WJ^ -VVij - — (£,,-+A,,--nu;)^ + (r2+7)^' 

n 

where state \i) and \j) are in the left and right well respec- 
tively, and £ij has the same definition as Eq. (4). Furthermore, 
Tqi = exp{ — (£12 — Ai2)/T}. When the frequency is 5 
MHz and the dephasing rate is 27r x 1 GHz, we use the or- 
dinary cooling method to obtain the population distribution 
in state |1) at different amplitudes and detunings as shown in 
Fig. 10 (a). Fig. 10 (b) is the contour plot of the population 
in state 1 1) with the new method in the same condition as that 
in Fig. 10 (a). Extracting the lowest population in Fig. 10 (b) 
and its corresponding amplitude at each detuning as shown in 
Fig. 10 (c) and Fig. 10 (d) respectively, we find the cooling 
efficiency in the new method is much better than that in the or- 




Microwave Frequency (MHz) Microwave Rux Amplitude (mOg) 

FIG. 1 1 : (Color online), (a) The population in state 1 1) versus the fre- 
quency and flux amplitude with the detuning 6^dc = 0.05 m$o and 
the dephasing rate r2 /27r = 1 GHz in the new method. The other 
parameters used are the same with those of Fig. 3. (b) and (c). The 
lowest population in state |1) extracted from (a) at each frequency 
and the corresponding amplitude versus microwave frequency. 



dinary method shown in Fig. 8 (b) especially at the detuning 
near m^o- Without the influence of A03, the amplitude con- 
dition for optimal cooling keeps consistency at all detunings 
and satisfies the condition 2<l>^/+ 6^dc = 8.646 m<l>o ^ ^12- 
Because the strong decoherence enlarges the transition region 
near A 12, the amplitude needs to be increased to longer the 
time in the region. 

In order to compare two methods, we put them together in 
Fig. 10 (a). The white solid line represents the amplitude 
for optimal cooling in the ordinary method, while the green 
dashed line represents that in the new method. Actually, the 
green dashed line is part of the optimal region in the ordinary 
method with different cooling interests shown by the white 
dashed arrowheads. Cooling interests of the green solid line 
have large energy spacing and hence the influence of the ther- 
mal noise is little. However the cooling interests in the new 
method are bothered by the temperature due to a small spac- 
ing. Therefore the green dashed line produces great effects to 
realize lower effective temperature on the cooling interests. 

At the detuning S^dc = 0.05 m<l>o we change the ampli- 
tude and frequency to obtain the population in state |1) with 
the new method, shown in Fig. 11 (a). Then by extracting 
the lowest population at each frequency in Fig. 11 (a), an un- 
conventional phenomenon occurs that larger frequency yields 
better cooling as shown in Fig. 11 (b). This phenomenon dis- 
obeys the frequency condition in the ordinary method, and we 
would discuss this result in the next section. In addition, in 
Fig. 11 (b), although 5 MHz is not the frequency for opti- 
mal cooling, it also provides an effective cooling and we can 
still use 5 MHz for cooling in practice. Moreover at each fre- 
quency the amplitude corresponding to the lowest population 
in state |1) as shown in Fig. 11 (c), demonstrates that larger 
amplitude is needed to produce the optimal cooling at larger 
frequency. Since larger frequency results into shorter time 
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FIG. 12: (Color online). The lowest population in state |1) among 
different amplitudes versus the dephasing rate and microwave fre- 
quency at the detuning 6^dc =0.05 m$o in the new method. The 
parameters used here are the same with those of Fig. 3. 

in the transition region, larger amplitude is needed to brings 
more time in the transition region over one period. However 
more transition region covered by the microwave indicates av- 
eragely smaller transition rate described by Eq. (8). These two 
factors result into the behavior in Fig. 11 (c). 

Furthermore, still at the detuning S^dc = 0.05 m<l>o, at 
different dephasing rates and frequencies we vary the ampli- 
tudes to obtain the lowest population in state 1 1), as illustrated 
in Fig. 12. The frequency corresponding to the optimal cool- 
ing increases as the dephasing rate increases instead of staying 
at 5 MHz. Since the new method eliminates the disturbance 
of Ao3 which conflicts with the cooling mechanism, the con- 
dition for optimal cooling here should agree with that of the 
experiment. This anti-intuitive result indicates there are un- 
derlying and interesting physical mechanisms. 

V. OPTIMAL MDLZ TRANSITION 

To explain the strange phenomenon in Sec. IV, firstly we 
return to MDLZ transition and change the transition rate Eq. 
(9) to the form 

TT/ _ (2u;M,,-)^/^Ai^[(2u;/A,,)^/^(n-A,,/u;)] 

^^^^ ~ 2 ^ (£^,+A,,■-nu;)V(^2+7)+^2+7 ' ^^^^ 
n 

where we use the airy function Ai(x) to express the Bessel 
function Jn{x) = (2/x)^/^Ai[(2/x)^/^(n - x)]. The sum is 
determined by the terms with n, the nearest integer to (Sij — 
Aij)/uj. Therefore, we have 

W.J « Al{^)iAi^[{2uj/A)h,,M/{T2 + 7)+ 

Ali^)iAi'[{2u;/A)He,, +u;)/co]/{^ + ^). 

(11) 

Then we would use the above expression to talce a quahtative 
analysis. Eq. (11) shows the term w^/^/(r2 + 7) + + 
7)/w^/^ would increase the second and third term until the 
frequency reaches to nearly (r2 + 7)/27r and in this process 
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the total transition rate would increase. This explains why 
under the dephasing rate T2/27T = 1 GHz the increase of fre- 
quency from 5 MHz to 500 MHz brings the decrease of pop- 
ulation in state |1) in Fig. 11 (b). As the frequency is larger 
than +7) / 27r, this term would contribute to the decrease of 
the transition rate Wij . Therefore the frequency correspond- 
ing to the maximum of Wij is relative with the decoherence 
rate. Since here the contribution from other term in Eq. (10) 
is neglected, we only give a qualitative analysis and further 
conclusions could not be obtained from Eq. (11). 

Therefore, we return to the expression in Eq. (9). For a 
small dephasing rate, e. g. r2/27r = 0.05 GHz, the frequency 
corresponding to the maximum of Wij would be nearly 10 
MHz [Fig. 13 (a)]. If the dephasing rate r2 is large, e. g. 
r2/27r = 1 GHz, the maximum of Wij would need a large fre- 
quency nearly 470 MHz, as shown in Fig. 13 (b). That means 
in this case the frequency of the maximum of Wij reaches at 
the incoherent transition region introduced in Sec. II, not stay- 
ing at nearly 5 MHz. Fig. 13 (d) is the plot of the frequency 
corresponding to the maximum of Wij vs. the dephasing rate. 
The frequency increases with the dephasing rate, agreeing 
with the qualitative analysis. MDLZ transition rate Wij de- 
pends on the driving period and the population transferred by 
LZ transition over one period. The two factors are both related 
with the microwave frequency. Larger frequency yields more 
times of LZ transitions in one unit time, but also brings less 
LZ transition population. The above qualitative and quantitive 
discussions demonstrate at an appropriate frequency the max- 
imum of MDLZ transition rate can be obtained. Furthermore 
the decoherence has an important effect on the modulation of 
the appropriate frequency. 

When the dephasing rate is small, e. g. r2/27r = 0.05 
GHz, at the detuning S^dc = 0.05 m^o the population in 
state |1) shown in Fig. 13 (c) does not at first decrease to a 
lowest value as W12 indicates in Fig. 13 (a), but exhibits a 
monotonic increase. This is because the population in state 
|1) near the degeneracy point is dependent on not only W12, 
but also 1^10 and W23. In this case the low frequency driving 
would result the low velocity near Aqi and A23 such that Woi 
and W23 have contribution on the population in state If 
we choose the detuning far away from the degeneracy point, 
this difference between the population in state |1) and W12 
would disappear as shown in Figs. 14 (a) and 14 (b). 

When the dephasing rate is large, e. g. r2 /27r = 1 GHz, the 
population exhibits a monotonic decrease with the frequency 
and does not increase again [Fig. 11 (b)]. At the detuning far 
away from m^o, the population in state 1 1) is 
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FIG. 13: (a) and (b). The maximum of W12 among different ampli- 
tudes at each frequency versus microwave frequency at the detuning 
S^dc = 0.05 m$o under the dephasing rate r2/27r = 0.05 GHz 
and 1 GHz respectively, (c) shows the lowest population in state 
|1) among different amplitudes at each frequency and the detuning 
S^dc = 0.05 m$o under the dephasing rate r2/27r = 0.05 GHz. 
(d). The frequency of the maximum of W12 versus the dephasing 
rate at the detuning 5^dc = 0.05 m$o- All figures are calculated 
in the new method. The other parameters are the same with those of 
Fig. 3. 



where Fq^ = F'^q exp{ — (£12 — Ai2)/T}, and Wij is defined 
by Eq. (9). Seen from Eq. (12), Fq^ would also play an im- 
portant role in the cooling. In order to produce the optimal 
cooling, larger frequency needs larger amplitude [Fig. 1 1 (c) 
and Fig. 15]. According to the expression of Fq^ in Eq. (12), 
larger amplitude would result the decrease of Fq^, which di- 
minishes the population in state Therefore, the final pop- 



ulation exhibits different characters from W12. 

To summarize, the decoherence modulates the frequency of 
the maximum of MDLZ transition rate, influencing the cool- 
ing effect. Although the incoherent transition region does not 
receive much attention in quantum information due to the loss 
of coherence, our studies show that the region has the par- 
ticular advantage in quantum transition and quantum cool- 
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FIG. 14: (a). The maximum of W12 among different amplitudes 
at each frequency versus the frequency in the new method with the 
detuning S^dc = 3 m$o and the dephasing rate r2/27r = 0.05 GHz. 
(b). The lowest population in state |1) at each frequency versus the 
frequency with the detuning 6^dc — 3 m$o and the dephasing rate 
r2/27r = 0.05 GHz. The other parameters are the same with those 
of Fig. 3. 
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FIG. 15: (Color online) The amplitude corresponding to the max- 
imum of W12 at each frequency versus the frequency in the new 
method with the detuning S^dc = 0.05 m<^o at r2/27r = 0.05 GHz 
(the solid line) and 1 GHz (the dotted line). The other parameters 
used here is identical with those of Fig. 3. The inset shows the am- 
plitude corresponding to the maximum of W12 at each dephasing rate 
and frequency with the detuning 6^dc — 0.05 m$o. 



ing. This decoherence induced modulation of the maximum of 
MDLZ transition should also emerge in the ordinary cooling 
method as indicated in Eq. (11). Under the weak decoher- 
ence, e. g. r2/27r = 0.06 GHz, at a detuning far away from 
Aoi, e. g. 5^dc = 0.5 m^o, even when the state is driven to 
completely reach the transition region near A 12, A03 would 
not have effects. Therefore, nonmonotonic effect of MDLZ 
transition can be found in the population just as shown in Fig. 
6 (b) with a frequency ~ 10 MHz corresponding to the opti- 



mal cooling. However, at the detuning 5^dc = 0.05 m^o the 
amplitude of the maximum of 1^12 makes the crossover A03 
effective. (see the inset of Fig. 15) Therefore, smaller ampli- 
tude is needed and results of a final frequency 5 MHz. 



VI. CONCLUSION 

Starting from a recent experiment, we investigate the 
microwave-driven cooling in an artificial atom subjected to 
the environment noise. Under the weak decoherence, we show 
that the optimal cooling at each detuning depends on not only 
the microwave frequency but also the amplitude. At each de- 
tuning, the frequency where the optimal cooling is achieved 
is not always at 5 MHz, maybe slightly larger than 5 MHz. 
At the same time, the optimal cooling also requires that the 
amplitude satisfies the relation ^rf ^ *^i2 — ^^dc- 

We further discuss the cooling under the strong decoher- 
ence, where more vibrational degrees of freedom are coupled 
to the two lowest states, making the ordinary cooling method 
ineffective. To recover the effective cooling, we present an 
improved cooling method, which employs a revised driving 
waveform to avoid another "harming" side crossover. The 
new method can be used to effectively improve the cooling ef- 
ficiency, prepare the quantum state and suppress decoherence 
in multi-qubit system under not only the weak decoherence 
but also the strong decoherence. This robustness of effective 
cooling under the decoherence is useful when the qubit is em- 
ployed to cool neighboring quantum systems and refrigerate 
the environmental degrees of freedom. In addition, small mi- 
crowave frequency needed in the new method also provides 
the possibility to use our method on adiabatic quantum com- 
putation. 

Furthermore, we show that better cooling effect is not de- 
pendent on lower frequency and the frequency for optimal 
cooling in both methods monotonically increases with the de- 
coherence stronger. Currently, MDLZ transition in the co- 
herent region has attracted lots of attentions due to coherent 
phenomena. We point out the incoherent transition region 
has the particular advantage in the cooling and MDLZ transi- 
tion with both methods. Our study demonstrates qualitatively 
and quantitively the frequency where MDLZ transition rate 
reaches maximum is modulated by the decoherence and deep- 
ens the understanding to the physical mechanism of MDLZ 
transition as well as the optimal cooling. Moreover, the re- 
sults based on MDLZ transition can be used to other systems 
with similar energy states structure, realizing quantum state 
manipulation. Especially, for systems consisting of the in- 
teraction between the qubit and the physical resonator, such 
as micromechanical beam, our study can further be used to 
achieve the optimal cooling of the qubit, or even resonators. 



VII. ACKNOWLEDGEMENTS 

We thank Valenzuela, S. O. for useful discussions. This 
work was partly supported by the State Key Program for 
Basic Researches of China (2011CB922104), the NSFC 



11 



(10725415), and the Natural Science Foundation of Jiangsu Province (BK2010012).. 



Electronic address: yuyang@nju.edu.cn 
^ Y. Makhlin, J. Diggins, and A. Shnirman, Rev. Mod. Phys. 73, 
357 (2001). 

^ M. A. Nielsen and I. L. Chuang, Quantum Computation and 
Quantum Information, (Cambridge University Press, Cambridge, 
2000). 

^ J. Q. You and F. Nori, Phys. Today 58, No. 1 1, 42 (2005); 

J. E. Mooij, Science 307,1210 (2005) 
^ J. Clarke and F K. Wilhelm Nature(London) 453, 1031(2008) 
^ Y. Nakamura, Y. A. Pashkin, and J. S. Tsai, Nature(London) 

398,786(1999) 

^ D. Vion, A. Aassime, A. Cottet, P. Joyez, H. Pothier, C. Urbina, 
D. Esteve, and M. H. Devoret. Science 296 886 (2002) 

^ Y. Yu, S. Y. Han, X. Chu, S. I. Chu, Z. Wang. Science 296, 
889(2002) 

^ J. M. Martinis, S. Nam, J. Aumentado, and C. Urbina. Phys. Rev. 
Lett. 89,117901(2002) 

I. Chiorescu, Y. Nakamura, C. J. P. M. Harmans, and J. E. Mooij. 
Science 299, 1869(2003) 

Y. Nakamura, Y. A. Pashkin, and J. S. Tsai, Phys. Rev. Lett. 87, 
246601 (2001) . 

I. Chiorescu, P. Bertet, K. Semba, Y. Nakamura, C. J. P. M. Har- 
mans, and J. E. Mooij, Nature (London) 431, 159 (2004) 
A. Wallraff, D. I. Schuster, A. Blais, L. Frunzio, R.-S. Huang, 
J. Majer, S. Kumar, S. M. Girvin, and R. J. Schoelkopf, Nature 
(London) 431, 162 (2004) . 

D. I. Schuster, A. A. Houck, J. A. Schreier, A. Wallraff, J. M. 
Gambetta, A. Blais, L. Frunzio, J. Majer, B. Johnson, M. H. De- 
voret, S. M. Girvin, and R. J. Schoelkopf. Nature (London) 445, 
515 (2007) 

A. Izmalkov, M. Grajcar, E. Il'ichev, T. Wagner, H.-G. Meyer, A. 
Yu. Smirnov, M. H. S. Amin, A. M. van den Brink, and A. M. 
Zagoskin, Phys. Rev. Lett. 93, 037003 (2004) 
M. Grajcar, A. Izmalkov, S. H. W. van der Ploeg, S. Linzen, T. 
Plecenik, T. Wagner, U. Hubner, E. Il'ichev, H.-G. Meyer, A. Y. 
Smirnov, P. J. Love, A. M. van den Brink, M. H. S. Amin, S. 
Uchaikin, and A. M. Zagoskin, Phys. Rev. Lett. 96, 047006 (2006) 
W. D. Oliver, Y. Yu, J. C. Lee, K. K. Berggren, L. S. Levitov, and 
T. P Orlando. Science 310, 1653(2005) 

M. Sillanpaa, T. Lehtinen, A. Paila, Y. Makhlin, and P. Hakonen. 
Phys. Rev. Lett. 96,187002(2006) 

A. Izmalkov, S. H. W. van der Ploeg, S. N. Shevchenko, M. Gra- 
jcar, E. Il'ichev, U. Hubner, A. N. Omelyanchouk, and H. G. 
Meyer, Phys. Rev. Lett. 101, 017003 (2008). 
C. M. Wilson, T. Duty, F. Persson, M. Sandberg, G. Johansson, 
and P Delsing, Phys. Rev. Lett. 98, 257003 (2007) 
J. Hauss, A. Fedorov, C. Hutter, A. Shnirman, and G. S. Schon, 
Phys. Rev. Lett. 100, 037003 (2008). 

M. Grajcar, S. H. W. van der Ploeg, A. Izmalkov, E. Il'ichev, H.- 
G. Meyer, A. Fedorov, A. Shnirman, and G. S. Schon, Nat. Phys. 
4, 612 (2008). 

I. Martin, A. Shnirman, L. Tian, and P Zoller, Phys. Rev. B 69, 
125339 (2004) 

M. Grajcar, S. Ashhab, J. R. Johansson, and F. Nori, Phys. Rev. B 
78, 035406 (2008) 

J. E. Mooij, T. P Orlando, L. Levitov, L. Tian, C. H. van der Wal, 
and S. Lloyd, Science 285, 1036 (1999). 

I. Chiorescu, Y. Nakamura, C. J. P. M. Harmans, and J. E. Mooij, 



Science 299, 1869 (2003). 

S. O. Valenzuela, W. D. Oliver, D. M. Berns, K. K. Berggren, L. 
S. Levitov, T. P Orlando, Science 314, 1589 (2006) 
D. J. Wineland, R. E. Drullinger, F L. Walls, Phys. Rev.Lett. 40, 
1639 (1978). 

W. Neuhauser, M. Hohenstatt, P. Toschek, H. Dehmelt, Phys. Rev. 
Lett. 41, 233 (1978). 

I. Marzoh, J. I. Cirac, R. Blatt, P Zoller, Phys. Rev. A. 49, 2771 

(1994) . 

C. Monroe, D. M. Meekhof, B. E. King, S. R. Jefferts, W. M. 
Itano, D. J. Wineland and P Gould, Phys. Rev. Lett. 75, 4011 

(1995) . 

S.N. Shevchenko, S. Ashhab, and F. Nori, Physics Reports 492, 
1 (2010). 

A. V. Shytov, D. A. Ivanov, and M. V. Feigel'man, Eur. Phys. J. B 
36, 263 (2003). 

L. D. Landau., Phys. Z. Sowjetunion 2, 46 (1932); G. Zener, Proc. 
R. Soc. London. Ser. A 137, 696 (1932); E. C. G. Stuckelberg, 
Helv. Phys. Acta 5, 369 (1932). 

J. Ankerhold and H. Grabert, Phys. Rev. Lett. 91, 016803 (2003). 
G. Ithier, E. Collin, P Joyez, D. Vion, D. Esteve, J. Ankerhold, 
and H. Grabert, Phys. Rev. Lett. 94, 057004 (2005). 
K. Saito, M. Wubs, S. Kohler, P Hanggi, and Y. Kayanuma, Eu- 
rophys. Lett. 76, 22 (2006). 

K. Saito and Y. Kayanuma, Phys. Rev. B 70, 201304(R) (2004). 
G. Z. Sun, X. D. Wen, B. Mao, J. Chen, Y. Yu, P H. Wu, and S. 
Y. Han, Nat. Commun. 1:51, doi: 10.1 03 8/ncomms 1050 (2010). 
S. N. Shevchenko, S. H. W. van der Ploeg, M. Grajcar, E. Il'ichev, 
A. N. Omelyanchouk, and H.-G. Meyer, Phys. Rev. B 78, 174527 
(2008). 

D. M. Berns, M. S. Rudner, S. O. Valenzuela, K. K. Berggren, W. 
D. Oliver, L. S. Levitov, and T. P. Orlando, Nature(London) 455, 
51(2008) 

G. Z. Sun, X. D. Wen, Y. W. Wang, S. H. Cong, J. Chen, L. Kang, 
W. W. Xu, Y. Yu, S. Y. Han, and P H. Wu. Apphed Physics Letters 
94,102502(2009) 

J. Bylander, M. S. Rudner, A. V. Shytov, S. O. Valenzuela, D. M. 
Berns, K. K. Berggren, L. S. Levitov, and W. D. Oliver, Phys. Rev. 
B 80, 220506 (2009) 

W. D. Oliver, S. O. Valenzuela, Quant. Info. Proc. 8, 261 (2009). 
D. M. Berns, W. D. Oliver, S. O. Valenzuela, A. V. Shytov, K. K. 
Berggren, L. S. Levitov, and T. P. Orlando, Phys. Rev. Lett. 97 
,150520(2006) 

L. J. Du, M. J. Wang, and Y. Yu, Phys. Rev. B 82, 045128 (2010) 

L. J. Du, and Y. Yu, Phys. Rev. B 82, 144524 (2010) 

M. S. Rudner, A.V Shytov, L. S. Levitov, D. M. Berns, W. D. 

Oliver, S. O. Valenzuela, and T. P. Orlando, Phys. Rev. Lett. 

101,190502(2008) 

A. O. Niskanen, Y. Nakamura, and J. P. Pekola, Phys. Rev. B 76, 
174523 (2007) 

J. Q. You, Y. X. Liu, and F Nori, Phys. Rev. Lett. 100, 047001 
(2008) 

M. H. S. Amin, and D. V. Averin, Phys. Rev. Lett. 100, 
197001(2008) 

^2 M. H. S. Amin and F Brito, Phys. Rev. B 80, 214302 (2009) 
L. J. Du, and Y. Yu, in preparation. 



